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Abstract 

In this paper, proceeding from the recently developed way of deriving the 
quantum-mechanical equations from the classical ones, the complete system 
of hydrodynamical equations, including the quantum Euler equation, is derived 
for a perfect fluid and an imperfect fluid with pairwise interaction between the 
particles. For the Bose-Einstein condensate of the latter one the Bogolyubov 
spectrum of elementary excitations is easily reproduced in the acoustic approx- 
imation. 
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Introduction 

Scientific interest in the de Broglie-Bohm causal interpretation of quantum 
mechanics [2, Q and its applications has been recently appreciably rekindled 
(see also, for example, In particular, in following the ideas of N.G. 

Chetaev (the theorem on stable trajectories in dynamics |5() and G. 't Hooft 
(the classical deterministic theory, supplemented with the mechanism of dis- 
sipation, generates the observed quantum behaviour of our world [6]), it was 
explicitly demonstrated, that adding the dissipation energy Q to the Hamilton 
function and imposing the stability condition on motion of the corresponding 
mechanical system, one can come to the standard Schrodinger equation, as well 
as identify the famous Bohmian potential exactly with the introduced dissipa- 
tion energy. Besides, the indisputable advantage of the produced derivation of 
the Schrodinger equation from Newtonian mechanics lies in the fact, that \ip\ 2 
turns out to represent the density of trajectories in the configuration space. 
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One of the next natural steps lies in developing the hydrodynamical approach 
to the problem of motion and deriving the quantum Eulcr equation, and then 
applying it for describing the Bose-Einstein condensates. This paper is devoted 
exactly to this very important and promising line of investigation. Obviously, 
if fundamental properties of quantum fluids are theoretically reproduced in the 
Bohmian mechanics in a natural way without any additional assumptions and 
ungrounded extensions, than this fact may serve as one more circumstantial 
evidence of this theory. 

The paper is organized in the following way. In Section 1 we derive the 
modified Liouville equation, taking into consideration quantum effects. Then 
we come to the complete system of the quantum hydrodynamics equations in 
Section 2. Finally, in Section 3 we consider sound waves in a quantum im- 
perfect fluid with pairwise interaction between the particles and reproduce the 
Bogolyubov spectrum of elementary excitations for its Bose-Einstein conden- 
sate. The main results are summarized in Conclusion. 



1. Modified Liouville equation in de Broglie-Bohm theory 

One of the most logical and self-consistent methods for deriving the macro- 
scopic hydrodynamical equations consists in averaging of the microscopic me- 
chanical equations of motion of each single particle. Following this prevalent in 
fluid physics method, let us consider a system of N identical particles, whose 
motion obeys Newton's second law: 



dri dvi dU 

— =Vj, m— = Fi(T k ,t), Fi(r fe ,i) = - — 
at at ori 



(1) 



where (t) and Vj (t) are the radius- vector and the velocity of the i-th particle 
respectively, m is the mass of a single particle, is the force, acting on the «-th 
particle from the direction of the external field (for example, the gravitational 
one) and all other particles, and U(vk,t) is the corresponding potential energy. 
For simplicity we restrict ourselves in this paper to the case of uncharged parti- 
cles in the absent external electromagnetic fielcfE Here and in what follows the 
index k corresponds to the set of natural numbers from 1 to TV inclusive. 

In order to proceed from the microscopic equations ([1]) to the macroscopic 
ones, one needs the TV-particle distribution function /^(r^, v/., t), satisfying the 
standard normalization requirement 

N 

/ /i\r(rfc, Vfc,t) JJdr^v; = 1 (2) 



lr The more general case of charged particles in the present external electromagnetic field is 
a subject of our forthcoming paper. 
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and the corresponding well-known continuity equation 
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in the phase space. Along with the function /jv(rfe, t), representing the 
probability density in the phase space, should satisfy the following evident equa- 
tion: 



N 

/ /iv(rfe,v fe ,i) JJdvi = A 2 (r fc ,i), 



(4) 



where the function A(rfc,i) = |-0(rfe,t)| is the modulus (amplitude) of the N- 
particle (total) wave function -0(rfc,i), obeying the Schrodinger equation 
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Obviously, as it directly follows from ^ and (U]), for the functions A(rk,t) 
and ip(jk,t) the standard normalization requirement 



j A 2 (r fc) t) Y[dn = I \^(r k ,t)\ 2 l[dr q = 1 



(6) 



holds true (naturally, A 2 (i*k,i) — \tp(rk,t)\ 2 represents the probability density 
in the configuration space). Let us also note, that the functions iW(rfc, Vfc, i) 
and A(rk,t) are invariant relative to permutation of two arbitrary particles in 
view of their identity. 

According to the de Broglie-Bohm theory [H, 0, H, Q , the potential energy 
U (rfc , t) may be presented in the form 



U(r k ,t) = U cl (r k ,t) + Q(r k ,t) , 



(7) 



where U c i(rk,t) is the classical part, and Q{vk,t) is the Bohmian quantum po 
tcntial: 

Expressing A(r^ , t) from Q and substituting the result into © , we obtain 



(8) 
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Substituting ((T|), ([7]) and © into (O, we come to the modified Liouville 
equation 
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Obviously, the modification lies in the last term of quantum nature, pro- 
portional to ti 2 . In contrast to three previous terms, it is essentially nonlinear 
concerning fN(r k , t). It should be also mentioned, that the obtained equa- 
tion (1101) contains all information about the system. It other words, strictly 
speaking, it is not necessary to solve any other additional equation, for exam- 
ple, the Schrodinger equation ([S]), for solving (fit)]), where all quantum features 
are already taken into account. 

On the other hand, the modified Liouville equation (fTO)) is very difficult to 
solve, so one needs to carry out some simplifications. The clearest one of them 
is considered in the next section. 



2. Equations of Bohmian hydrodynamics for a perfect fluid 

Let us give concrete expression to the classical potential energy U c \(r k ,t), ne- 
glecting pairwise interaction between the particles (or, in other words, regarding 
our "fluid" as a perfect one): 

N 

U cl (r k ,t) = Y,V(r j ,t), (11) 

i=i 

where V{tj , t) is the potential energy of the j-th particle in the external field. For 
simplicity let us assume, that all N particles are " in the same quantum state" , 
described by the normalized to unity wave function with the modulus Aq, 
then the solution of the Schrodinger equation ([3]) may be presented in the form 

iir k ,t) = l[Mr J ,t), ih^-=(~A + V(r,t)\^o, (12) 

2=1 

whence it follows, in particular, that 

N 

A(r k ,t) = l[Ao(r j ,t). (13) 

2=1 
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It should be mentioned, that the function Aq(t, t) represents the macroscopic 
density of an arbitrary single particle. Consequently, the macroscopic density 
n(r,t) of the whole system is interconnected with Aq(t, t) as follows: 



n(r,t) =NAl(r,t) A (r,t) = 

Substituting (JT3J) and (0 into ©, we get 

h 2 A, Agfa,*) 



n(r, t) 



N 



(14) 
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(15) 

Let us note, that the made assumption of the same quantum state for all 
N particles enables to present the iV-particle quantum potential Q(rfc,t) in the 
form of the sum of N identical 1-particle ones Qi(rj,t), each depending on the 
radius-vector Yj of the corresponding j-th particle, j = 1, 2, . . . , N. Simultane- 
ously the iV-particle amplitude A{r^,t) is presented in the form of the product 
of N identical 1-particle ones Ao(rj,t). A similar presentation holds true for 
the A-particle distribution function: 



/iv(rfc,Vfc,t) = Y[fi(rj,Vj,t) 

3=1 



(16) 



where the 1-particle distribution function /i(r, v, t) satisfies the normalization 
condition 



J fi(r,v,t)drdv 



= 1 



(17) 



As it directly follows from (|10p after the substitution of (fTTj) and (|16p (or 
from (j3)) after the substitution of (JTJ, (jlip and (|15p). this function obeys 
the equation 
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Applying the standard methods of fluid physics (in particular, presenting 
v in the form of the sum of the hydrodynamical velocity field u(r, t) and the 
deviation <5v from it), one can easily obtain from (fl"8|) the complete system of 
hydrodynamical equations, describing motion of our perfect fluid: 
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(19) 
(20) 
(21) 
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where the mass density p(r, t), the pressure p(r, t), the temperature T(r, t) and, 
finally, the heat flux density q(r, t) are the macroscopic quantities, given by the 
following equations respectively: 

p(r, t) = mn(r, t) = Nm J fc (r, v, t)dv , (22) 

p(r, t) = n(r, t)reT(r, t) = ^iVm / Ju 2 A(r, v, t)dv , (23) 

q(r,i) = ^Nm J Sv 2 Svf 1 (r, v, t)dv . (24) 

Here k is the Boltzmann constant, and viscosity is neglected. The continuity 
equation (jTSJ) and the heat equation (f2~Tj) exactly coincide with the correspond- 
ing equations of classical (non-quantum) hydrodynamics Q . Only the equation 
P0|) , which represents the quantum Euler equation, differs from its classical 
analogue in the presence of the last term in the right hand side, proportional, in 
particular, to the ratio h 2 /m 2 . Besides, it is absolutely obvious, that one should 
expect the greater absolute value of the quantum potential Q, the sharper the 
spatial change of the mass density p is. If the fluid is incompressible, that is p 
is simply some constant, then Q = 0, and the distinction between classical and 
quantum Euler equations is eliminated. In the general case the additional quan- 
tum term in (|20[) vanishes, if the function J~p satisfies the Hclmholtz equation 
A yfp + const^/p = 0, where const is an arbitrary constant. 

Let us estimate, for which typical velocities and distances in some concrete 
problem this term, being nonzero, plays a vital part. Introducing typical velocity 
Mo and distance ro, one can immediately compare it with the second term in 
the left hand side of the same equation. They become quantities of the same 
order for UqTq ~ h/m. In this simple relationship one can easily recognize the 
famous Heisenberg's uncertainty principle, as it should be. For example, if it is 
a question of helium, then the product UqTq should be a quantity of the order 
1CT 8 m 2 /s. 

Thus, in this section we have derived the equations (fl"9]l . (|20)) and (|2ip 
of Bohmian hydrodynamics for a perfect fluid. Concluding this section, let us 
note, that in the framework of standard quantum mechanics the hydrodynamical 
equations were written down by Landau and Khalatnikov (see, for example, [7(), 
but they contain quantum-mechanical operators with proper commutation rules 
and are much more complicated. 



3. Sound waves in imperfect fluid with pairwise interaction 

According to [1,0, only for temperatures, close to absolute zero, quantum 
effects advance to the forefront in properties of fluids, but actually only helium 
remains liquid down to absolute zero, while all other fluids become solid before 
quantum. As regards helium, our basic assumption of the same quantum state 
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for all particles holds true for the Bose-Einstein condensate of its isotope 4 He, 
being a Bose fluid (see also [HI 12 1). Now we are interested in developing 



quantum acoustics for such quantum fluids, being, generally speaking, imperfect 
owing to interparticle interaction. 

In order to take proper account of pairwise interaction between the particles, 
one can add in the right hand side of (TTTI) the term 



1 N 

V ?i=2 E $ ( r ;- r '), (25) 

where $(i"j — r;) is the potential energy of interaction of j-th and l-th particles, 
and then apply the method of correlation functions for solving (JTofE but there 
is also another way. Even in spite of interparticle interaction, the form of the 
equations from (| 13[) to (|16[) may be preserved as a certain approximation, if the 
" 1-particle" Schrodinger equation (fT2"]) is modified in the following way: 

lh ^dT = (~£ A + y(r ' t] + VGp{r > t] ) ^ ' (26) 

where the additional potential energy 

V GP (r, t) = NA 2 (r, t) = —p(r, t) , (27) 

m m 

where, in its turn, a is a typical length of boson-boson scattering. The nonlin- 
ear 1-particle Schrodinger equation ([26)) is well-known as the Gross-Pitaevskii 
equation 12, 13| (see also 13, HI)- Taking into consideration the presence of 



27|) . instead of ([20)) we get the quantum Euler equation in the form 

du ( d\ I dp 1 dV 47rh 2 adp h 2 d ( ' Ay/p\ 

— +u—u = —/--- JL + i^LL\ 28 

at \ or J p or m or ra A or 2m z Or \ Jp J 

At the same time two other hydrodynamical equations (11^1) and (|2T|l remain 
unchanged. 

Now let us turn to the case V = and apply in the equations (flU]) and 
(|28[) the standard acoustic approximation: up to the zero order of smallness 
p = po, u = 0, p = pq (where po and po are some constants, depending neither 
on time t, nor on spatial coordinates r), and up to the first order of smallness 
p = pa + pi, u = ui, p = po + pi (where the additional terms pi, ui and p\ 
are some functions of t and r, as well as the presence of the adiabatic relation 
Pi = c 2 Pi is assumed, where c s is the sound speed). Then from (TIT?)) and (|2"8"j) 
up to the first order of smallness we obtain respectively 

dpi dui , . 



2 Development of this method for the modified Liouville equation in the de Broglie-Bohm 
theory is also a subject of our forthcoming paper. 
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<9ui c 2 dpx 4irh 2 a dp 1 h 2 d 
ot po or m i Or Am z pa or 
From (|29|) and (|3U|) the sound propagation equation follows immediately: 

d 2 Pl ( 2 47rfi 2 op \ /i 2 



<9f 2 \ s m 3 ) Am 2 

Let us establish the dispersion law for plane waves. For that let us look for 
the solution of (l3Tj) in the form 

pi(r, t) ~ exp[i(kr) — iujt] , (32) 

where k is a wave vector, and a; is a frequency. Substituting (|32[) into (1311) . we 

where = |k|. It should be mentioned, that here the sound speed squared c 2 
is proportional to temperature and, consequently, becomes negligible for tem- 
peratures, close to absolute zero, while the scattering length a depends on tem- 
perature weakly, remaining nonzero. Then the established dispersion law (|33l) 
exactly coincides with the famous Bogolyubov spectrum of elementary exci- 
tations in a Bose fluid. Thus, this fundamental spectrum is derived in our 
quantum-hydrodynamical approach in a sufficiently simple way. 



Conclusion 

In this paper, proceeding from the recently developed way of deriving the 
quantum- mechanical equations from the classical ones [4j , we have derived suc- 
cessively the modified Liouville equation (10) and the complete system of hy- 
drodynamical equations (namely, the continuity equation (19), the quantum 
Euler equation (20) and the heat equation (21)) for a perfect fluid, formed by 
identical particles, being in the same quantum state. This basic assumption 
holds true for Bose-Einstein condensates of quantum fluids at very low temper- 
atures. Then we have taken into account the interparticle interaction by using 
the Gross-Pitaevskii equation (26). Finally, using the acoustic approximation, 
we have easily reproduced the Bogolyubov spectrum of elementary excitations 
(see (33)). ' " _ _ 

In future, being based on the direct relativistic generalization [14| of [J], 
we shall also generalize the developed here quantum hydrodynamics to the rel- 
ativistic case. This line of investigation promises to open up new intriguing 
possibilities for studying, in particular, black holes and the whole Universe. 
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